function [] = fig_sec10_2_adaptiveRoutingDual()

close all;

xp = load ('data/adaptiveRoutingDual_xp.txt');
pie = load ('data/adaptiveRoutingDual_pie.txt');
ye = load ('data/adaptiveRoutingDual_ye.txt');
objFunc = load ('data/adaptiveRoutingDual_objFunc.txt');
jomReg_xp = load ('data/adaptiveRoutingDual_jomReg_xp.txt');
jomReg_pie = load ('data/adaptiveRoutingDual_jomReg_pie.txt');
jomReg_ye = load ('data/adaptiveRoutingDual_jomReg_ye.txt');
jomReg_objFunc = load ('data/adaptiveRoutingDual_jomReg_objFunc.txt');

multiplot ('Time', 'Regularized cost' , objFunc , jomReg_objFunc , [80 150] , '-' , '\pi_e' , pie , jomReg_pie , [0 1.1] , '-' , 'y_e' , ye , jomReg_ye , [0 30] , '-' , 'x_p' , xp , jomReg_xp , [0 30] , '-');

% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-900 -500] , '-' , 'p_e' , pe , jom_pe , [0 0.2] , '-' , 'u_e' , ue , jom_ue , [0 0.1] , '-');
% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-100 100] , '-' , 'p_e' , pe , jom_pe , [0 1] , '-' , 'u_e' , ue , jom_ue , [0 1] , '-');

jomReg_objFunc
grad_objFunc = objFunc (end)
